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Abstract 

We compute annihilation amplitudes for charmless B decays that are proportional to the three- 
parton twist-3 light meson distribution amplitude 4>3m{xi, $2) with an active gluon. Due to an 
enhancement from a quark propagator at the scale p 2 ~ to;,Aqcd these terms occur at the same 
parametric order in a s (mb) and 1/mb as the known leading order annihilation involving fs and 
twist-2 meson distributions. With our calculation the leading order annihilation amplitude is now 
complete. At lowest order in a s the new amplitudes are real and only O5-8 contribute. Using 
simple models we find that the three-parton and two-parton terms are of comparable size. 
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The nonleptonic charmless decay channels B — > M\Mi provide a wealth of information 
about the standard model, including the study of CP violation and the strong interactions. 
Since many amplitudes for these decays are loop dominated, it is possible for new physics to 
give a significant contribution. However, except for the simplest observables, testing for new 
physics requires an understanding of the standard model background. Predicting standard 
model decay rates and CP asymmetries with quantum chromodynamics (QCD) is difficult, 
but the task is simplified by the use of the soft-collinear effective theory (SCET) Q, Q, S, 
and factorization theorems 

Baai. 

An interesting experimental observable is the relative "strong" phase between standard 
model amplitudes multiplying the CKM factors V u bV*j and V c bV*j (/ = d,s), since these 
phases are measured to be large in the B — » nn and B — > Kn channels 0]. There are two 
competing standard model explanations for these phases, sizeable charm penguin loops 0, 
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12] or sizeable annihilation amplitudes 0, Q, 0, 0] in which the initial state 
"spectator" quark is Wick-contracted with a quark field in the effective Hamiltonian. In 
this paper we report on a new contribution to the leading annihilation amplitudes. 

Our notation follows that of Ref. 0] where factorization theorems for the leading order 
B — > M\M.2 amplitudes were derived with an expansion in A/Q where A is a hadronic scale 
and Q ~ m& ~ m c ~ Em- We restrict our discussion to non-isosinglet mesons (M, — tt, K, 
p, . . .) which can not be produced solely by gluons, for which the annihilation amplitudes 
are power suppressed by ~ A/Q. Recently these power corrections were classified according 
to their perturbative order and source for strong phases [171]. The annihilation ampltude is 

Aann _ Alann) , Alarm) A2arin) /, \ 

71 ~~ ^hard T ^hard-col T ^hard "+-•••> I 1 J 

where the superscript denotes the order in A/Q. A subscript "hard" denotes annihilation 
amplitudes generated by propagators offshell by p 2 ~ m 2 . At lowest order in a s these 
include the standard leading order amplitudes A^^ 1 ~ [«s(^6)/b/j\/ 1 /a/20a/i0a/ 2 ] as well 
as the chirally enhanced amplitudes . \, 2 'Z'' : ~ ot-s (mh) fn Im, Im^ 4> m-^ <$m 2 to,/ J . These 



have been studied in earlier analyses Ja [13|, [14|, [15|, and at this order, factorization in 
rapidity 18J reveals that they are real [17j. Here the /'s are decay constants, and the chiral 
enhancement factors are fi n = m 2 T /{m u + m^) and fix = m 2 K/{ m u + m s)- </>Mi are twist- 2 
distribution functions, and <^ is a two-parton twist-3 distribution. In Eq. ([I]) the amplitudes 
^hard-coi nave hard-collinear propagators, which are offshell by an amount jj 2 ~ m b A. A 
non-perturbative phase is first generated by soft exchange between the two mesons as an 
order a 2 (jj,i)/ir suppressed term in Aj^_ col . The ellipsis in Eq. (OQ) denotes the fact that 
the full set of ^( 2ann ) amplitudes have not yet been classified. 

In this paper w( 
which has the form 



In this paper we compute the leading term in the perturbative expansion of A^^j^, 



^hard-col ~ a s {m b ) H ^ Xl, ^ 1,y ^ ® f B (f>B( k ) Sm^Mj. («l) /sA/ 2 03Af 2 (j/l, V2) ■ (2) 
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Here if is a calculable hard-scattering kernel, 03m is a three-parton twist-3 distribution, and 
f 3 M is the corresponding decay constant. The amplitude in Eq. <^ occurs at the same order 
in I / mb and a; s (m&) as v4.[ i ^ ra ' ) and should be included for a complete leading order anni- 
hilation amplitude. Unlike A^™^ its convolution integrals converge without using rapidity 
factorization. Furthermore, the LO annihilation involves 5-meson information beyond fg, 
thus demonstrating that annihilation is more complicated than the short distance picture 
leading to a scaling ~ fB/ m b that is often used in parametric estimates [l9j . 

In QCD at the scale m&, flavor changes are mediated by the weak effective Hamiltonian. 
For B -> M X M 2 with AS = 0, 



rr _G F 



10,77,89 

E fa°? + C *°2 + E c *°*) ■ ( 3 ) 



p=u,c 



Most of these operators have spin (V—A) ® (V— A), such as 0" = (w6)v_^ (dit)v-yi. We will 
prove below that all such operators give vanishing contribution to Eq. (j2j), so that only 

°5 = Y. q '{db) v ^A (qq)v+A , 6 = Y. q <{dp h a )v-A (qW/sh+A , 

°7 = Eg'-^- (WW , 8 = Y, q >-£- (dpb a ) v „ A {q' a q'p) v+A , (4) 

are relevant for our analysis. Here a and (3 are color indices, e q r are electric charges, and 
the sum is over flavors q' = u, d, s, c, b. Results for AS = 1 transitions are obtained by 
replacing d — > s in Eqs. ((31) and (jlj), and likewise in the equations below. The coefficients in 



Eq. ([3]) are known at next-to- leading- log order 20] . (We have 0\ <-> Of relative to [20]). In 
the NDR scheme with m b = 4.8 GeV, the coefficients are C 5 _ 8 (m fe ) = {0.010, -0.040, 4.9 x 
10~ 4 , 4.6 x 10~ 4 }, They are considerably smaller than Ci(mj,) = 1.08 and ^(mj) = —0.18, 
but can give important contributions in penguin observables because only contribute 
through loops [141 ] . 

To separate the mass scales below mj we match Hw onto operators in SCET. The am- 
plitude for B — > MiM 2 is most easily calculated in the B rest frame where soft fields with 
typical momenta ~ A interpolate for the initial state B. The final state hadrons Mi and M 2 
are back-to-back energetic charmless pseudoscalar or vector mesons. Collinear fields in the 
light-like direction n interpolate for one light meson, and collinear fields in the direction n 
interpolate for the other. These fields have typical momenta (n-p,n-p,p±) ~ Q{ff-, 1, rj) and 
Q(l,r] 2 ,r]), respectively, in terms of the power counting parameter rj ~ A/Q. The vectors 
n and n satisfy n ■ n = 2, and we work in a frame where the S-meson four- velocity v has 
n-v = n-v = 1. To calculate the amplitude in Eq. §2§ we first match QCD onto operators in 
SCETj where the hard-collinear modes with p 2 ~ m^A are still propagating and then match 



these onto operators in SCETn which has only non-perturbative modes with p 2 ~ A 2 [21]. 

Before presenting the details of the calculation of Eq. we complete our review of 
Ref. [l3] with a discussion of the required ingredients and an overview of the matching 
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FIG. 1: Picture for Q^C^, which generates an annihilation amplitude that is sensitive to the 
intermediate scale p 2 ~ m^A. The filled circle at the center represents a SCETi six-quark operator 
arising from the full-theory diagrams in Fig. (|2|) at the hard scale p 2 ~ m 2 . 

procedure. Let denote the leading SCETi weak operators. Annihilation contributions 
require either (a) an n-quark n-antiquark pair produced by hard interactions giving a 

Q(k>2) 

six-quark operator, or (b) a time-ordered product of two mixed Lagrangians C^ q C^J q to 
produce the n-n pair by soft exchange. Annihilation also requires a mechanism for connecting 
the "spectator" to the weak operator, either by (i) having a soft field directly in the weak 
operator, or (ii) having a time-ordered product with an Cf^ . Case (a,i) gives operators 
which contribute only to A^^ . Case (b,i) vanishes at this order. Case (b,ii) involves three 
hard-collinear gluons and is ~ a 2 (/j,i). This leaves case (a,ii). Here Q^C^ q can contribute 
at 0(a s (m b )) if the gluon from C} q is uncontracted. Since the uncontracted gluon costs an 
extra power when matched onto SCETn, it is only this class of operators that contributes 
with an external gluon, not Q^C^ q , Q^C^ q \ etc. Thus the calculation of A^rd-coi involves 
finding SCETi operators of the form 

where Q us <S> 6n <S> O n are color and spin structures, q and q' are flavors, and the collinear 
direction n' = n or n. The fermion fields are gauge invariant with large label momenta 
specified by the subscripts uj, for example q-n,^ = 5(uji — n-V)WlC,n ' ) where W n is a Wilson 
line. At tree level these operators arise from the full-theory diagrams in Fig. [2] with three 
light n'-collinear quarks and two collinear in the other direction, n'. They have Wilson 
coefficients of 0(a s (rrib))- We identify n as the collinear direction of the pair-produced 
quark of flavor q and sum over all n in the SCETi weak Hamiltonian. We will see shortly 
that the flavor structure is as in Eq. and that the matching requires T A color structures 
for two of the 0's. 

To finalize our description of the calculation we consider matching the time-ordered 
product T[Q^Cf q ] ] onto SCET n with dia grams as shown in Figure [H has an ex- 

cess of n'-collinear fermions since only two are needed to interpolate for a collinear me- 
son. The sub leading Lagrangian [22] C^ q = q' us ig$n'l'n> removes an n'-collinear fermion 
and provides the soft field that interpolates for the light anti-quark in the B meson. Here 
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FIG. 2: Tree-level annihilation graphs for B — ► M1M2 decays. The gluon and the fermion propa- 
gator connecting it to the weak vertex are both offshell by p 2 ~ m&. Matching on to SCETj, these 
graphs give rise to the six-quark operators Q^ 2 \ the filled circle at the center of Fig. [TJ 

igB^ = [\/{n'-V) Wl[in ■ D n ,,iD», ± }W n ,5(uj - n' and the form of the SCET n 
operators is 

OfP oc — — [q s ^ n , k T s b v ] [d n T n q n ] [q n T n q' n ]igB^ , (6) 

with T s £g> T n ® T n containing spin and color structures. The collinear gluon field strength 
igB^, ~ rj, interpolates for gluons in a final state meson, so there is no perturbative suppres- 
sion from the factor of g. At tree level, integrating out the hard-collinear quark propagator 
in Fig. [1] induces an inverse factor 1/ (n' ■ k) of the soft momentum which will be convoluted 
with the ^-distribution, ^{n' -k). In Eq. ([6]) this compensates the 77 suppression from igB^> 
to make 0^ 1T ^ the same order as the six-quark operators for the hard annihilation, which is 
0{rf). We have checked that operators with more igB^'s or with soft gluon field strengths 
do not occur at this order in 1/rrib and a s (m(,). 

Note that SCETn time-ordered products (T-products) do not contribute at 0(r/ 7 ). To 
see this, recall that our process has a soft initial state and n and n-collinear final states. An 
example of an SCETn Lagrangian that connects these sectors 23( has two-collinear quarks 
and two-soft quarks [24], q s q s in£,n ~ V- 111 these operators the two n-collinear particles 
conserve the large p~ ~ rj momenta, and the two soft particles conserve the p + ~ r] 
momenta. Thus this operator, as well as analogous operators with gluons, only support 
scattering, ns — > ns, and not annhilation such as nn — > ss or ss — > nn. Another example is 
£11 ~ ^nAnAfi^n ~ rj 2 , where analogous statements hold for n and n. Weak operators, like 
0^ T \ that have the same n-n-s structure as the initial and final states are already 0(r) 7 ), 
so T-products with them are power suppressed. The above considerations rule out the 
majority of T-products. An example of an annihilation T-product in SCETn that survives 
these criteria is Cn, with a weak operator with fields (q s h v ^n^,n) ~ V 5 - These T-products 
involve at least one loop momentum where, due to the double multipole expansion, (r 
must be smaller than the conserved p~ and p + , see Eq.(25) of Ref. [2^. As a contour integral 
in £ + or £~ we have > 2 poles that are all on the same side of the axis, and therefore the 
loop gives zero. At 0(rf) this is sufficient to rule out possible annihilation T-products, 



including those with more than one SCETn Lagrangian. Note that in Ref. [24J a T-product 
contribution was identified for B° — > D°tt°, however in that scattering process the integral 
did not satisfy the same pole criteria as we find here. 
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Constructing the Operator Bases 

Next we construct a full basis for the operators and 0^ lT ^ in the SCETi and SCETn 
weak effective Hamiltonians, respectively. (The matching calculations beginning on page 7 
can be understood without the details of this somewhat technical construction, the results of 
which are Eqs. (Tl5!) and (fl6l) .) General symmetry arguments allow us to reduce the operator 
bases to the small subset relevant to our calculation of ^hMd-coi' anc ^ f° r this reason it is 
convenient to construct the bases for SCETi and SCETn simultaneously. First consider 
spin in SCETi. For light fermion fields of definite handedness, a complete basis of Dirac 
structures for the individual bilinears in Eq. §5$) is 

Qus/w = {i,72}, &n = {i,til}, e n = {f£,^}. (7) 

Using these bases, we must construct a complete set of spin structures with chiralities 
inherited in perturbative matching from the full-theory fields in Oi_io and the produced qq 
pair. To make a Lorentz scalar, the spin structure must have zero 7j_'s or two 7j_'s contracted 
with g°f = g a/3 — n a n^ /2 — n l3 n a /2. Note that contracting with e°f = n p n a e af3p,J /2 does 
not yield an independent operator since for example i^^n^lv^n = £n$7±75£n = £n$7x£n • 
For Oi-^g^o the only allowed chiral structure is (LH)(LL)(LL) where L and R refer to the 
handedness for the light quarks in the bilinears in the order shown in Eq. ([5]). We cannot 
assign a handedness to the heavy quark denoted here by H. This chiral structure is realized 
as the spin structures 

©us/n ® @n ® ©n = 1 ® f ® ft, ©us/n ® ©n ® @n = 1 ® f ® # • (8) 

We have ruled out the chirality (LH)(LR)(RL) corresponding to a spin structure 1 ®ifi^± ® 
0-7a by using P R ^ 7" ® Pi$la — 0. This equation encodes the helicity flip argument 
of Ref. [15j. Similarly, for Os_8 the chirality (LH)(RR)(RR) is also realized as the spin 
structures Eq. (|S]), whereas (LH)(RL)(LR) is not allowed since -Pl$7" <8> Pr^7^ = 0. We 
will show momentarily, however, that using SCETn the terms in Eq. (jSJ) are not needed to 
compute ^ard-coi- F° r ^5-8 we can also have 

<3>us/n <S> ©n ® ©n = 7± ® J* ® fl7a , ©W« ® ©n ® ©n = 7± ® ?^7a ® $ , (9) 

corresponding to chiralities (RH)(LL)(LR) and (RH)(LR)(RR), respectively, and thus the 
flavor structure shown in Eq. ([5]), namely (q'b)(dq)(qq'). The second structure in Eq. (jSJ) 
is related to the second structure in Eq. (|2J) by a Fierz transformation swapping d n and 
q' n quarks and we will choose the latter for our operator basis. The complete set of spin 
structures in Eqs. (jSJ) and Qj contains neither Q us / n <8> On ® ©n = 1±_ ® fila ® $ nor 
©us/n ® @n <E> O n = 7± ® ^ ® ^7a • These possibilities are excluded by the projection relation 
©us/n' = $> ^'©us/n'^/4 and the helicity flip equation. 

Now consider spin and chirality in SCETn- The allowed 0^ lT " > spin structures must 
respect the handedness inherited from the SCETi fields in the perturbative matching of 
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T[Q^Cf>}. For n' = n in Q^ 2 \ taking either one of the n-collinear anti-quark fields soft 
yields an annihilation operator. For n! = n, however, the field q n in the third bilinear was pair 
produced and does not contribute to the annihilation amplitude when made soft by C^ q . So 
given the SCETi spin structures Eqs. (IE]) and (J9j) corresponding to chiralities described in the 
text, we need to consider chiralities (LH)(LL)(LL), (LH)(RR)(RR), (RH)(LL)(LR), 
and (RH)(LR)(RR) with bilinears in the order shown in Eq. ([H]), i.e. soft — n — n. With the 
first bilinear purely soft, a complete basis of Dirac structures for the individual bilinears is 

r s = {^,7?}, r* = {i,ii3, r n = {^,^}. (10) 

A Lorentz scalar 0^ 1T ^ has an odd number of 7j_'s since one must be contracted into the 
n— or n— collinear field strength fij. For chiralities (LH)(LL)(LL) and (LH)(RR)(RR) the 
allowed Dirac structure is 

(r s ® T n <g> T n )B^ ± = (^®i/i® $)B p n , ]X (11) 

with n' = n or n, but the corresponding operators 0^ 1T ' have qsJ^K and do not contribute for 
B decays. Since (LH)(LL)(LL) is the only O^ 17 '^ chirality corresponding to the (V— A)(V— A) 
operators Oi_4 j9) io, this proves that only 5 _ 8 can contribute to Eq. (|2J). Furthermore since 
all (LH) terms are ruled out, the soft quark can only be q', and not a <i-quark. 

This leaves the (RH)(LL)(LR) and (RH)(LR)(RR) structures from Os_g with soft quark 
flavor q', for which we have the additional spin structures, 

ri = n : T s ® T n ® T n i3f ± = ® ifi ® #^ n ,±, ® ^ ® #i3f >± } , 

n' = n: ^(8)^5^®^ = {fl®^n^L®^ f^®^,±®^}, (12) 

plus those with ^ <-> fl in r s . While these eight are all allowed by chirality and Lorentz 
invariance, six can be ruled out by considering the spin and factorization properties of our 
time-ordered product. The matching from SCETj to SCETn does not affect the spin and 
color structure of the n'-collinear bilinear at this order in the power expansion, since once 
a jet direction is chosen the collinear fields in the opposite direction are decoupled. Here 
ft! is the opposite of n'. From Eqs. (JSD and dHD the allowed Q n i structures have no 7_i_'s, 
and therefore the second structure on each line of Eq. ffl2l) does not appear at any order 
in the perturbative matching. Also, the allowed structures Eqs. (jSj) and Q are invariant 
under Q us US $ /2 and only power-suppressed interactions couple the b— quark to the n' 
sector. Therefore, T s should not vanish under T s r s / /2, and the operators with iji «-> fi 
mentioned below Eq. (fT2l are ruled out. In perturbation theory this just corresponds to the 
appearance of an ifi from the n'-collinear propagator next to the 6-quark. This leaves only 
the operators with a.B± in Eq. ( TT2"j) . 

Finally consider color. In SCETi the operators are color singlets, but each bilinear 
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on its own could be singlet or octet. A complete set of color structures includes 

®us/n> ® ©n ® ©n = {T* ® 1 ® T°, T° ® T° <g) 1, 1 ® 1 ® 1, 

1 <g> T a £g> T a , T a ® T 6 <g> T c / a6c , T a <g> T 6 <g> T c d abc } . (13) 

Once again we can reduce this set using the factorization properties of SCETi. As argued 
for spin, an SCETi operator with color structure 0^ matches onto a SCETn operator with 
the same structure T n i in its n' bilinear. So Q n > cannot be a color octet, and the allowed 
structures are 

Qus/n ® ©« ® ©n = {1 ® 1 ® 1, T a eg) 1 <g> T a } 

©« s /n ® ©« ® ©„ = {1 ® 1 ® 1, T a ® T a ® 1} . (14) 

In SCETn each of the three bilinears interpolates for a color singlet meson and therefore 
each bilinear must seperately be a color singlet, F s eg) eg T n — 1 eg 1 eg) 1. 

Matching onto SCETi and SCETn 

We now present the matching from Hw in Eq. ([3]) onto the SCETi operators and 
then the matching of the SCETi time-ordered product T[Q^C^ q ] onto SCETn operators 
0^ T \ The hadronic matrix elements of wm gi ye the factorization formula for A^^_ col . 
From the arguments presented above, the complete basis of SCETi operators is 

m b 

Q?d — ~3 J2 q ,q' [^PLllT a K] [d n ^T a P R q n ^} [q n ^P R q' n ^] , 
m b 

(2) _ (2) 3ey 

with sums over q, q = u, d, s, plus analogous operators Q 5d _ 8d which have color structure 

(2) 

1 Cg) 1 Cg) 1. The electroweak penguin operators 07,8 induce the two operators Q 3d4d , which 

(2) 

have the same spin and flavor structures as 0\J 2d , but with a factor of the quark electric 
charge e q i included under the summation. Combining the pieces in SCETn, a complete basis 
for the 0{rf) operators with one igB± that contribute to B decays is 

1 
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°iP = — 3T+ H,-k+ P LiS ] n b v ] [d n ^P L q njUJ3 ] [q n ^$ (ig$ ± ) n:UJ5 P R q' nuJ4 ] 
°2d ] = -^rrEov [€,-h- p Li Slb v ] [dfl lWa # (ig$±)fi,u s PRqn,u 3 ] [qn^PRq'n^] 



mi k~ ^ q ' q ' 
(iT) _ (iT) 3e^ 

w 3-4d — w l-2d 2 ' ^ ' 

Here <^ = (q' s S n )8(k + + n-V') and q' s _ k - = (q' s Sn)S(k~ +n-V') and the direction for the 
soft Wilson lines S n and S n are determined by the matching from SCETi. Just like the local 
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(IT) 

annihilation operators, we see that the 0\ s can not create transversly polarized vector 
mesons. The basis for AS = 1 decays, switches d n — >• s fl . 

Next, we carry out the perturbative matching onto the bases in Eqs. (fl5l) and ffl6l) . and 
derive the factorization theorem. The SCETi weak Hamiltonian with Wilson coefficients af c 

(2) 

for the operators Q id is 

Hw = ^ (Af + Af ) £ / [cfc^dwadws^dws] £ a? fa) Qfjfa) . (17) 

* n,n i=l—8 

Since only the penguin operators Os_g contribute, we pulled out the common CKM factor 
with Au = VubVud an d A^ = V c bV* d . The analogous result for AS = 1 has the same a^ c 
coefficients. To match onto the af c at tree level we first do a spin Fierz on the full theory 
O^s operators to obtain spin structures Pl <8> Pr, and then compute the graphs in Fig. [21 
Only graphs c) and d) are nonzero and we find [at fi = mj 

h 7ra s (m b ) \ 2C F C 5 + C 6 (2C F - C A )C 5 + C 6 

ai ^ mV) N c \y[x(l-y)-l] + (l-x)y(l-y) 
n hc (T = ,_ na s (m b ) f (2C F -C A )C 5 + C 6 2C F C 5 + C 6 \ 

21 ' ' V) ~ N C \ S[(l-*)(l-y)-l] xy(l-x) ^ 



The coefficients a^ c 4 are identical to a\ c 2 respectively with the replacements C^^ — > C> 



7,8- 



-he 
*5-8 



also begin at 0(a s {mb)) but give «<; (/^-suppressed contributions when matched onto 
SCETn, so we do not list their values. These coefficients are "polluted" in that one- loop 
0(a s (mb) 2 ) contributions proportional to Ci 2 could compete numerically with the results 
in Eq. (fl8|) . Here x,x,y, and y are defined in Fig. [2j namely y = ui/rrib, y = —u^/nib, 
x = uj 2 /mb, x = —u 3 /mb. For n! = n as in 0^3, we have x = 1 — x, but y ^ 1 — y since the 
momentum is shared between three n— collinear partons. Likewise, for n! = n as in a 2 4 we 
have y — 1 — y but x ^ 1 — x. 

Having constructed the operators and determined their Wilson coefficients, it is 
straightforward to match the time-ordered products T[Q^C^} onto the SCET n operators 
0] T . For odd indices i and even indices i' we find that integrating out the hard-collinear 
quark propagator, shown as the dashed line inside the gray region in Fig. [TJ gives 

i J*xT[Q<j${ Uj )]{0)£VW ~ J dk+oSPik+w), 

1 J#xT[Q® fam^x) ~ J dk-0SP(k-, Uj ). (19) 

At 0(a 2 s ) in perturbation theory this matching would include non-trivial jet functions. For 
example, in the first line a J du' l A J(k + , uj' 14 , 001,4) with taking the place of ui^ in of d . 
However at this order additional time-ordered products and non-perturbative functions be- 
come relevant so we stick to 0(a s ) in our analysis. Together Eqs. (I18|19l) complete the 
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(IT) 

tree-level matching. Now take the matrix element of O id using 

I j j 

(Kn^lun^fiPi^Rdn^lO) = 5 nm S(n -p- U X +UJ 4 ) (f>p(y) , (20) 

(pnAPi^n.uAPLfidn^) = HYl^Ll. g S(fi- p-^+U^ (f>vAy) , 

Zn-p 11 

and the three-body distributions 

if 

(^n 1 (p)\Un,oj 1 ^(ig$±)n,u J5 PRdn,u J4 \0) = <W 5(fi -p -U X - U 5 + ^4) 3P (y, y) , (21) 

(p+ ,{p,e)\ w„ )£Jl $ (ig$±) n ,u J5 PRd n:UJ4 10) = l ^ vmvn § 5{n-p-u l -u 5 +u A ) <p 3 v(v,v) ■ 

co>5 n ■ p 

Our convention for the vector meson matrix element has been chosen to simplify the final 
result for the amplitude and is related to that of 26] by / 3 y = m v fy and 3 y = — T/2. 
Permutations in the flavors give the definitions for other meson channels, and we use the 
phase convention in j^]. The soft matrix element is 

ml%. k PLi' Sl,b v \B) = i f -^d>%{n' -k) . (22) 
Combining these pieces the factorization theorem with tree-level jet functions is 

i4£52 ~' = < A « ,+A ?'>i* nr (23) 



1 r l 
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x S faMifMa dx dy dy — — — </>Wi(v, y)4>M 2 (x 



i-v-y 

i ,i H MlM2 (x,x,y) 
+ VmJmJsm 2 j dy j dx j dx hc ^_ — — (p Ml (y)hM 2 (x, x) 



Jo Jo 



where rjM = —1 or +1 for a pseudoscalar or vector meson, respectively. The hard coeffi- 
cients Hf? c l m and Hj^ m for different B -> M X M 2 channels are listed in Table [I] in terms 
of coefficients in the SCETi weak Hamiltonian. The amplitude contains the three-body 
distribution function as promised. The convolutions in Eq. (I23I) are real, and assuming the 
standard endpoint behavior for the distribution functions they converge without the rapidity 
factorization of jl8| . 

We conclude by comparing our result parametrically and numerically to A^^j and 
^hard as defined in Ref. 17]. For this comparison it is useful to define moment parameters 

& 2 = Jdxdydy MvrfMx) , (24) 

Ptf\ & M2 = [dydxdx tuMMx, x), (3 B = U ftUk) , 
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B -► MiM 2 


Hhci 


Hhc2 


n°7T~ , /9 7r~ 7T p _ 

7T~7r , p _ 7r° 7r~p° 
7T + 7r~, Tr + p~, p + ir~ 

7T ^ , A° 

^o^(*)O j ^(*)0^0 


1 r ,hc 1 n hc 
73 °1 73 °3 
1 „/lC 1 1 ,,/ic 

7f a i +73 a 3 

r.hc 1 1 ~/lC 

—a x -t- 

„/ic 1 „hc 
a l ~ 2 a 3 

_fec l_/ic 
"l ~~ 2 a 3 

of + a| c 


1 „/lC 1 1 _/ic 

73 a 2 +73 a 4 

1 „hc 1 „/ic 

~73 a2 73 ° 4 

a 2 2 a 4 

-of + i of 
-of + |af 
-af - af 


7T K\ ' : p K u 

7T A w j p 1\ 

ir°K^ ,p°K 


of + af 

1 „/ic 1 „/ic 

"73 ai 73 a 3 

1 „hc 1 „hc 

V2 1 2V2 3 

_ n hc 1 1 n hc 
u l ~ 2 "3 


(2-2 ^4 

73 a 2 +73 a 4 
~73 a 2 + 373 4 

2 2 4 


5 S -» MiM 2 


Hhcl 


Hhc2 


A f , A 7T , J\ p 
K°7T°, K*°7T°, K°p° 


„hc i 1 „/ic 
-O x + 

1 „ftc 1 „/ic 

73 a i 173 a 3 


„hc l k 
2 2 4 

1 „hc j_ 1 „/ic 

7f a 2 + 373 a4 


K+K~, K*+K-, K+K*- 
K°K°, K*°K° , K°K*° 


„ fee 1 1 „fec 

— a x -1- 

_fec 1 _hc 
"l ~~ 2 a 3 


he 1 n hc 
a 2 2 a 4 

-af + i a f 



TABLE I: Hard functions for the annihilation amplitude Aj, in Eq. (|23p for B°, and -B s 
decays. The result for B~ — > 7r°7r~ is obtained by adding the results using the entries from the 
first two rows, and so vanishes in the isospin limit. 



where (3b = A^ 1 /3 has mass dimension —1. First we compare the leading-power annihilation 
amplitudes in B —>■ n + K~ . Dropping terms proportional to the tiny Wilson coefficients CV_8j 
we have 

R , +Ir) = AtZU- + K-) ^(Ag+Ag) ^[-fMiZ + hKf^m 

Parametrically, the moments in R\ have (3^ ~ (3hd ~ 0(a s (mb)), and the power counting 
of the prefactor is J^k^bI Ik ~ 1- Also there is no suppression from the hierarchy in the 
Cj's since f3± u involves C3, and C3 ~ C5 ~ Cq. Thus, we have shown that for consistency in 
the a s and l/m& expansion, the contributions A^^} col need to be included with the local 
contributions in the leading annihilation amplitude. Similarly we can compare the 

new hard-collinear annihilation amplitude to the chirally enhanced annihilation contribution 
in B° — > n + K~ . Isolating the terms proportional to the large coefficients C5 and Cq we have 
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Parametrically /3 X ~ flhc ~ a s W, and R(n + K~) ~ mB/3g/ 3l / '(fn/J-n) ^ rrib / fi n nib / A as 
expected. 

We conclude with a brief numerical analysis of the ratios R\ and R 2 . The Cj's are 
quoted below Eq. and we use a s (m b ) = 0.22, f K = 0.16 GeV, and f n = 0.13 GeV. 
fs = 0.22 GeV is taken from a recent lattice determination [28], the three-body decay 



constants f S K — 4.5 x 10 3 GeV 2 and / 37r ~ 4.5 x 10 3 GeV 2 come from QCD sum rules [29 ] 



and Pb — l/(.4GeV) was determined in a fit to nonleptonic data [30J. To model the 



nonperturbative meson distributions we truncate the conformal partial wave expansions 31] 



as 



<P M ( X ) = 6x(l -x)[l + af (6z - 3) + 6af (1 - 5a; + 5a; 2 )] , 
(f) 3M (x, x) = 360xx(l - x - xf [I + ^{7(1 - x - x) - 3}] . (27) 

Eq. ( I24p has convergent convolution integrals for these distribution functions. To estimate 
the moments f3 and the ratios R we vary the coefficients in Eq. (1271) in a conservative range 
inferred from recent lattice results |32| for the af f 's and QCD sum rules [29( for the u>3m's. 
Specifically we take aj = 0, af = 0.05 ± 0.02, a^ K = 0.2 ± 0.2, and w 37T)K = -1 ± 1. A 
Gaussian scan of the model parameters gives 

= -1.4 ± 0.4, /?f c f = 0.3 ± 0.1, R^n+K") = 0.3 to 1.2 , 
^ = -1.4 ±0.5, ^ = 0.1 ±0.1, R 2 (ir + K-) = -0.1 to 0.1 . (28) 

The denominators of Eq. (|25[) and fT26|) can vanish, giving large departures from Gaussian 
statistics. So for Ri and R 2 we quote the range that contains an equivalent number of 
points as one standard deviation for a Gaussian distribution. Eq. ( 128]) demonstrates that 
numerically the three-parton contributions to y 4( lann ) could be of the same size or larger than 
the local piece A^^\ Numerically, vtibPb]^! {j-K^i^) ~ 0.2 causing some suppression in 
R 2 {^ ~K~). It would be interesting to examine the size of these three-parton contributions 
in the fc^-approach of Ref. 0]. 

In this paper we computed the final missing term of the leading order annihilation am- 
plitude in B — > M\M 2 decays. These terms involve a three-parton distribution and need to 
be included for a complete analysis of annihilation. 



Acknowledgments 

We thank Z.Ligeti for discussion. This work was supported in part by the Office of Nuclear 
Physics of the U.S. Department of Energy under the cooperative research agreement DOE- 
FC02-94ER40818 (C.A. and I.S.), and DOE-ER-40682-143 and DEAC02-6CH03000 (I.R.). 



12 



I.S. was also supported in part by the DOE OJI program and by the Sloan Foundation. 



C. W. Bauer, S. Fleming, and M. E. Luke, Phys. Rev. D63, 014006 (2001), [hep^ph/0005275[ 
C. W. Bauer, S. Fleming, D. Pirjol, and I. W. Stewart, Phys. Rev. D63, 114020 (2001), 
|hep-ph/0011336| 

C. W. Bauer and I. W. Stewart, Phys. Lett. B516, 134 (2001), |h^ph70107001 



C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D65, 054022 (2002), [hep^ph/0109045. 
M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda, Nucl. Phys. B591, 313 (2000), 



hep-ph/0006124 



Y.-Y. Keum, H.-n. Li, and A. I. Sanda, Phys. Lett. B504, 6 (2001), |hep-ph/000 4004. 



J.-g. Chay and C. Kim, Phys. Rev. D68, 071502 (2003), hep-ph/0301055. 

C. W. Bauer, D. Pirjol, I. Z. Rothstein, and I. W. Stewart, Phys. Rev. D70, 054015 (2004), 



hep-ph/0401188. 



J. Charles et al., Eur. Phys. J. C41, 1 (2005), URL |http : //ckmf itter . in2p3 . fr7[ 



M. Ciuchini, E. Franco, G. Martinelli, and L. Silvestrini, Nucl. Phys. B501, 271 (1997), 
hep-ph/9703353} 

P. Colangelo, G. Nardulli, N. Paver, and Riazuddin, Z. Phys. C45, 575 (1990). 
A. R. Williamson and J. Zupan (2006), |hep-ph/0601214| 



C.-D. Lu, K. Ukai, and M.-Z. Yang, Phys. Rev. D63, 074009 (2001), hep-ph/0004213 



M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda, Nucl. Phys. B606, 245 (2001), 
|hep-ph/01041~luj 

A. L. Kagan, Phys. Lett. B601, 151 (2004), [hep^ph/0405134[ 

A. Khodjamirian, T. Mannel, M. Melcher, and B. Melic, Phys. Rev. D72, 094012 (2005). 



C. M. Arnesen, Z. Ligeti, I. Z. Rothstein, and I. W. Stewart (2006), (hep^ph/0607001 , 

A. V. Manohar and I. W. Stewart (2006), |hep-ph/06 05001. 

B. Blok, M. Gronau, and J. L. Rosner, Phys. Rev. Lett. 78, 3999 (1997), [hep-ph/9701396 
G. Buchalla, A. J. Buras, and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 

C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D67, 071502 (2003), [hep^ph/0211069[ 
M. Beneke and T. Feldmann, Phys. Lett. B553, 267 (2003), [hep^ph/0211358[ 

R. J. Hill and M. Neubert, Nucl. Phys. B 657, 229 (2003), |hep-ph/02li018j 

S. Mantry, D. Pirjol, and I. W. Stewart, Phys. Rev. D 68, 114009 (2003), [hep^ph/0306254[ 

I. W. Stewart (2003), [hep^ph/0308185[ 

A. Hardmeier, E. Lunghi, D. Pirjol, and D. Wyler, Nucl. Phys. B682, 150 (2004). 

M. Gronau, O. F. Hernandez, D. London, and J. L. Rosner, Phys. Rev. D50, 4529 (1994). 



A. Gray et al., Phys. Rev. Lett. 95, 212001 (2005), (hepd^/0507015 . 



P. Ball, V. M. Braun, and A. Lenz, JHEP 05, 004 (2006), |hep^ph /0603063. 



C. W. Bauer, I. Z. Rothstein, and I. W. Stewart (2005), |hep-ph/051024l] 

A. R. Zhitnitsky, I. R. Zhitnitsky, and V. L. Chernyak, Sov. J. Nucl. Phys. 41, 284 (1985). 



V. M. Braun et al. (2006), hepdat/0606012 



13 



